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© Method for identifying complex lithologies of subsurface formations. 



© A method of identifying complex lithologies of 
subsurface rock formations includes the excitation of 
at least a portion of the rock formation with a multi- 
frequency electric current. The resulting voltage in 
the rock is measured and the phase and amplitude 
of such voltage are determined for each of the 
included frequencies. The resistivities are deter- 
mined for each of the frequencies and are plotted as 
a function of frequency. The complex lithology of the 
rock formation is characterized by such plot. 
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METHOD FOR IDENTIFYING COMPLEX LITHOLOGIES OF SUBSURFACE FORMATIONS 



The invention relates to a multiple frequency 
electric excitation method for identifying complex 
lithologies of subsurface formations. 

In the search for hydrocarbons and in the eval- 
uation of coal-and synfuei-type deposits, drilling 
represents an expensive commitment. A drilled 
hole can become worthless unless basic logs taken 
in such a hole can provide information vital to 
evaluating the potential of the surrounding area. If 
the lithology of the area is simple and known, 
analysis of the logs will give good values of the 
parameters needed to evaluate a hydrocarbon res- 
ervoir and to determine reserves. However, the 
complex composition of mixed lithologies being 
encountered worldwide in oil and gas exploration 
and production has greatly complicated log analy- 
sis. Also, current log analysis methods are not 
adequate in many areas. There is, therefore, a 
need for logs to specifically and quantitatively iden- 
tify basic rock types (limestone, dolomite, sand- 
stone, and mixtures of these), clay type and vol- 
ume, and minerals that adversely affect the logs 
used to estimate reserves. 

An electric log is a recording of the electrical 
resistivity of the rock formations surrounding a 
borehole. The measured resistivity, or electrical 
conduction, has been related to water saturation by 
the equations set forth in "The Electrical Resistivity 
Log As An Aid In Determining Some Reservoir 
Characteristics", Trans. AIME, Vol. 46, pp. 54-62, 
1942, by G. E. Archie. These equations, which 
have been so useful in the analysis of rock forma- 
tions, are often referred to as "Archie's Laws". 
Such equations are strictly applicable to the analy- 
sis of clean formations, that is formations with a 
matrix of framework that is so nearly an insulator 
that total electrical conductance is determined by 
conduction through the tortuous paths of the pore 
water. However, not many clean rocks exist in 
nature and log analysts recognized that a strict 
application of the Archie equations to shaly rock 
analysis produced a higher than actual water satu- 
ration. More recently, analysis methods have in- 
cluded electrical conductance of the clay minerals 
found in the shaly rock in a reformulation of Ar- 
chie's equations (see "Electrical Conductance In A 
Porous Medium", Geophysics . Vol. 48, No. 9, pp. 
1258-1268, Sept. 1983, by A. E. Bussian, and 
"Electrical Conductivities In Oil-Bearing Shaly 
Sands", Spa Petr. Eng. J., Trans. Vol. 243, pp. 
107-122, 1968, by W. H. Waxman and L J. M. 
Smits). 

However, clay minerals are not the only min- 
erals that can add to the total electrical conduc- 
tance of a rock. Electrically conductive minerals 



such as pyrite-and graphite-like organic matter can 
contribute and even dominate electrical conduc- 
tance in a rock. If this conductance is not properly 
taken into consideration, the log analyst can easily 

5 mistake a low resistivity and virtually water-free oil 
reservoir for one that is water saturated. The 
present invention seeks to provide a new method 
for measuring the electrical characteristics of rocks 
having conductive matrix elements, such as pyrite- 

w and other graphite-like organic matter, as well as 
clay minerals. 

The present invention is therefore direction to a 
method for identifying complex lithologies of a rock 
formation and, more particularly, to the identifica- 

75 tion of low resistivity rock formations that may be 
water-free oil or gas reservoirs. 

In accordance with that method, the rock for- 
mation is excited with a broad band multi-frequen- 
cy electric current. The resulting voltage in the rock 

20 formation is measured and the phase and am- 
plitude of such voltage at each of the included 
frequencies is determined. The resistivity of the 
rock formation is determined at each of the fre- 
quencies from the phase and amplitude determina- 

25 tions and is plotted as a function of such fre- 
quencies. The rock formation may then be char- 
acterized as (i) a clean sandstone when the resis- 
tance remains constant and the phase is negligible 
as a function of frequency, (ii) a low resistivity, high 

so surface area, clay-bearing rock when the real part 
of the impedance (resistance) remains constant 
and the exciting phase is a negative value which 
decreases in magnitude as the current frequency 
increases, and (iii) a low resistivity pyrite-bearing 

35 rock when the real part of the impedance - 
(resistance) decreases as the exciting current fre- 
quency increases and the phase is a negative 
value which increases in magnitude as the exciting 
current frequency increases. 

40 More specifically, the multi-frequency electric 

current excitation may be by means of white noise, 
by means of a plurality of discrete sine waves 
which are simultaneously applied, or by means of 
frequency modulation with a plurality of discrete 

45 sine waves. The frequencies utilized are suitably of 
a broad band width, for example from about 0.001 
hertz to 20 kilohertz. 

The invention is described below in greater 
detail by way of example only with reference to the 

so accompanying drawings, in which 

FIGS. 1-8 illustrate the resistivity and phase 
response of various complex lithologies to a broad 
band simultaneous frequency excitation; 
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FIGS. 9 and 10 illustrate logging systems for 
recording an electric log of rock formations sur- 
rounding a borehole; 

FIG. 1 1 illustrates a system for recording an 
electric log of a rock sample taken from a subsur- 
face formation; 

FIGS. 12-14 are detailed circuit schematics 
of three different embodiments of the current gen- 
erator of FIGS. 9 and 11 for generating the mul- 
tifrequency sine waveform used in the method of 
the invention; and 

FIGS. 15-17 are detailed circuit schematics 
of the current, voltage and phase detectors, re- 
spectively, of FIGS. 9 and 1 1 . 

In accordance with the method of the invention, 
rock resistivity is characterized in terms of its resis- 
tance and phase-angle as a function of a broad 
band frequency excitation of the rock. This is in 
contrast with the fixed, discrete frequency excita- 
tion currently in practice. The basis for a broad 
band simultaneous frequency excitation and a re- 
sponse in terms of rock resistance (real part of the 
impedance) and phase-angle (between excitation 
signal and response signal) as a function of such 
frequency excitation is illustrated in FIGS. 1-6 of 
the accompanying drawings for rock samples taken 
at different discrete frequencies over about five 
decades of frequencies ranging from about 0.01 
hertz to 200 hertz. These FIGURES clearly show 
that clay-bearing rocks and pyrite-bearing rocks 
have much different frequency-dependent complex 
resistivity responses. 

FIGS. 1 and 2 show that a clean rock, like a 
Berea sandstone, reacts to current excitation like a 
pure resistance. The phase-angle is very small, 
less than t 0.5 milliradian over the frequency 
range. FIGS. 3 and 4 show that a high surface 
area, clay-bearing rock (smectite) produces a mea- 
surable phase-angle that is negative and decreases 
in magnitude as the frequency increases. FIGS. 5 
and 6 show that a rock combining dispersed pyrite 
produces a measurable phase angle that is nega- 
tive and, in contrast with the clay-bearing rock of 
FIGS. 3 and 4, increases in magnitude as the 
frequency increases. In addition, this resistivity of 
this pyrite-bearing rock is frequency dependent; 
that is, the rock resistivity decreases as the excita- 
tion frequency increases. 

Additional FIGURES, FIGS. 7 and 8, depict a 
synthetic rock sample where two different types of 
electrical conduction takes place in parallel, one 
being purely electronic, such as continuous pyrite 
or graphite laminae, and the other being ionic con- 
duction by way of the pore water network in the 
rock sample. The rate of conduction is much faster 
in the purely electronic condution path as com- 
pared to the parallel ionic conduction path. The 
result is a large measurable positive phase-angle at 



the lower excitation frequencies which decreases 
as the higher frequencies are approached. There is 
also a frequency dependance in the measured 
resistance which suggests that a transient is being 
s measured rather than a true equilibrium type of 
phase shift. This synthetic rock represents a real 
case wherein a rock contains a conductive mineral, 
such as pyrite, in the form of an electrically con- 
tinuous lamina or coating on the walls of the pore 

io fluid network. 

These eight FIGURES are illustrative of the 
phase-angle and resistivity dependence on fre- 
quency in select rocks. Such dependency suggests 
that a subsurface formation surrounding a borehole 

75 can be uniquely characterized as to its' complex 
lithology by the running of the broad band multi- 
frequency electric log. This will eliminate water 
saturation measurement errors in clay-bearing and 
in pyrite-and other graphite-like-bearing rocks, 

20 which otherwise result if such frequency depen- 
dence is not taken into account, thereby resulting 
in the correct identification of low resistivity oil 
producing reservoirs which might otherwise be by- 
passed as water-saturated reservoirs. 

25 In carrying out the method of the invention, the 

rock whose complex resistivity is to be measured 
is electrically excited from a current source consist- 
ing of a broad band of frequencies. The measure- 
ment may be made in-situ by means of a borehole 

30 logging tool to produce an electric log of the rock 
formations surrounding the borehole or may be 
made on core samples taken from such rock for- 
mations. When the broad band multi-frequency 
current is injected into the rock, the voltage re- 

35 sponse of the rock is measured. This broad band 
multi-frequency current contains a plurality of fre- 
quencies in the range of near zero to about 20 
kilohertz, and consequently, excitation of the rock 
with such current is at all the frequencies simulta- 

40 neously. The amplitudes and phases of each fre- 
quency are known and all frequencies are prefer- 
ably in phase. The particular complex lithology of 
the rock transforms the current signal input into a 
unique voltage "fingerprint" consisting of dominant 

45 frequencies, amplitudes and phases which are dif- 
ferent from that of the input current. This voltage 
response and the transform that produces it 
uniquely identifies the rock in terms of its electrical 
properties. Another way of looking at the rock is as 

so a frequency-pass filter where certain frequencies 
are dominant in the measured response. 

A logging tool for carrying out cdmplex resistiv- 
ity measurements for identifying the complex lithol- 
ogy of subsurface rock formation, surrounding a 

55 borehole is illustrated in FIG. 9. The logging sonde 
10 employs a dipole-dipole electrode array. The 
broad band frequency excitation current from the 
current source 15 is induced into the rock forma- 
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tions through the pair of current electrodes 11 and 
12, and the resulting voltage is measured across 
the pair of voltage electrodes 13 and 14. These 
voltage electrodes are connected to a voltage de- 
tector 16 which measures the differential voltage 
across such electrodes. The differential voltage 
measurement is compared to the current from 
source 15 by detector 17 to determine the phase- 
angle. The voltage and phase are recorded in the 
form of the electric log 18. From these recorded 
measurements, the correct electrical resistance and 
reactance of the rock formation can be determined. 
An alternative arrangement, termed the "Wenner 
Array", for the current and voltage electrodes is 
illustrated in FIG. 10. While the dipole-dipole array 
of FIG. 9 will yield the best resolution, the Wenner 
Array will yield the best sensitivity. Other arrays 
are also possible, the arrays of FIGS. 9 and 10 
being merely two examples. More details as to the 
use of such a borehole logging tool with a single 
discrete frequency excitation current source for in- 
duced polarization logging may be obtained from 
U.S. Patent 4,464,930. 

A laboratory analysis on core samples taken 
from subsurface rock formations may be carried 
out with the system illustrated in FIG. 11. A core 
sample 20 is stimulated with a broad band fre- 
quency excitation current from the current source 



25. The voltage differential across the core sample 
20 is measured by the voltage detector 24 at 
terminals 22 and 23 which are located away from 
the ends of the cores ample so as to eliminate 

5 possible electrochemical errors resulting from the 
current input connection through resistor 28 and 
the current output connection to ground. Current 
through the core 20 is measured across a series 
resistor 28 by the current detector 29. The mea- 

10 sured voltage and current are compared by detec- 
tor 21 to determine the phase-angle between such 
voltage and current. This voltage current and 
phase-angle are recorded on the electric log 27. 
From these recorded measurements, the correct 

J5 resistance and reactance of the core sample can 
be determined. More details on the use of such a 
laboratory system for core sample electrical re- 
sistivity measurements are given below in conjunc- 
tion with the description of FIGS. 15-17. 

20 The band width of the broad band frequency 
excitation current required for carrying out the 
method of the invention will generally be in the 
range of about 0.001 hertz to about 20 kilohertz. 
The low frequency end is determined for a 

25 borehole logging operation by the logging speed 
and the depth increment of the measurement 
points. In general, the low frequency end of the 
band width is given by: 



f(hz) = LS (m/min) 



For example, assuming a logging speed (LS) of 12 solute phase, but also the measurement of resistiv- 

m/minute and a depth increment (A p) of 0.5 m, ity as a function of frequency. Illustrated in FIG. 12 

the low frequency end of the band width for 1 full are, for example, five sine wave generators 30-34 

cycle will be 0.4 hertz. producing five discrete signals (f, -f 5 ) at 1 Hz, 10 

The broad band frequency excitation current Hz, 100 Hz, 1 kilohertz and 10 kilohertz, respec- 

may be of three types: simultaneous discrete fre- tively. The outputs of each of these generators are 

quencies, frequency modulation, or white noise. applied through resistors 40-44 to resistor 45 and 

Referring firstly to simultaneous discrete frequency summing amplifier 36 to produce one waveform I- 

current excitation, such current excitation does not (t) that is a mixture of these five signals. If the sine 

include all frequencies in a selected band width, wave characteristic (sin at) of each discrete signal 
but rather, well-defined frequencies in the band 45 is the same, and the amplitudes (a) of each is the 

width. This type of excitation has many advantages same, the resulting signal l(t) takes the form: 
which include not only the measurement of ab- 



Kt) 



u 5 



sin (Jt + s0j) 



(2) 
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where A = gain x amplitudes and 
0 = phase-angle. 

The voltage response of the rock V(t) can be 5 



V(t) = 3 sin (Jt + 0j) . (3) 



Allowing for a DC component I, and V, in both l(t) 75 
and V(t), respectively, equations (2) and (3) can be 
rewritten as: 



I(t) = I Q Aj sin (Jt * i.), and (4) 



J=u, 



3j sin (Jt + 



(5) 



The l„ A J( Bj, V 0 , 9 U and 0 v values can be 
obtained from a matrix solution of N simultaneous 
equations, where N is the number of points of 
digitization of both the composite l(t) and V(t) 
waveforms. An analysis of V(t) gives both the am- 
plitude and phase at each of the discrete fre- 
quencies and logs can be plotted of phase versus 
frequency and resistivity versus frequency. 

Secondly, frequency modulation current excita- 
tion is similar to simultaneous discrete frequency 
current excitation and differs only in that there is a 
sweeping through the frequency band width rather 
than a simultaneous generation. Illustrated in FIG. 
13 is a circuit for carrying out such a frequency 
modulation excitation comprising a voltage ramp 
generator 130 and voltage controlled oscillator 131. 

Thirdly, white noise current excitation is the 
combination of all frequencies in the select band 
width at the same average amplitudes. The unique 
characteristic of the rock is determined by exciting 
the rock with white noise and measuring the fre- 
quency content of the induced rock voltage. Cer- 
tain of the input frequencies will still be present and 
certain frequencies will be absent. Those frequen- 
cies that still are present will be changed in phase 
and amplitude. Such an elimination of certain fre- 



quencies and the changing of the phase and am- 
plitude of the other frequencies identifies the trans- 

35 fer function that is characteristic of the rock. Illus- 
trated in FIG. 14 is a circuit for carrying out such a 
white noise excitation. A white noise diode 133, 
along with resistor 134 and capacitor 135 provide a 
wide flat spectra of noise. This noise is applied 

40 through amplifier 136 to band pass filter 137. The 
fitter 137 eliminates noise outside of the band of 
interest. 

The detailed operation of the current detector 
29, voltage detector 16 and 24, and phase detec- 
ts tors 17 and 21 of FIGS. 9 and 11 will now be 
described in more detail in conjunction with FIGS. 
15-17. In the current detector of FIG. 15, the volt- 
age across series resistor 28 of FIG. 11 is mea- 
sured by the amplifier combination 47 and 48, 
50 amplifier 48 functioning as a differential amplifier 
and amplifier 47 functioning as a buffer amplifier. 
These amplifiers are biased by resistors 49-52. The 
output of amplifier 48 is an analog representation of 
the excitation current flow through the rock sample. 
55 In the voltage detector of FIG. 16, the voltage 

differential from across the electrodes 13 and 14 of 
FIG. 9 or the terminals 22 and 23 of FIG. 11 is 
applied through preamplifiers 60 and 61 to voltage 
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differential amplifier 62. Amplifier 62, along with 
amplifiers 68-70, produce a desired voltage range 
for a given current input and rock impedance en- 
countered. Four ranges are illustrated for rocks of 
impedances up to 100 ohms, 1 kiiohm, 10 kilohms 
and 100 kilohms. The particular scale may be 
selected by way of switch 71 . Bias conductions are 
set by resistors 76-85, and the output of switch 71 
is an analog representation of the voltage of the 
rock formation or rock sample. 

In the phase detector of FIG. 17, the current I 
is applied to a fixed DC time delay provided by 
variable resistor 100 and capacitor 101 and then to 
amplifier 102. The voltage V is applied to the 



amplifiers 103 and 104. The output of amplifiers 
102 and 104 are summed to provide the phase 
comparison signal 0 on line 105. This signal is 
filtered through amplifier 106 to provide a phase- 

5 angle signal on line 107 representing the amount of 
phase shift in the measured voltage signal created 
by the reactive set by resistors 108-118, capacitors 
1 20-1 21 , and diodes 1 22-1 24. 

It should be understood, however, that the cir- 

10 cuitry described above with reference to FIGS. 9- 
17 is merely representative of numerous possible 
embodiments of the invention. The following are 
specific types of circuit components that may be 
used: 

75 



20 



25 



30 
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Reference Designation 
Oscillators 15, 25 and 30-34 
Oscillator 131 
Ramp Generator 130 
3and Pass Filter 137 
Amplifiers 47, 48, 60-62, 58-70 

and 136 
Amplifiers 36, 102 and 104 
Amplifiers 103 and 106 
Diodes 122-124 
Diode 133 

Capacitor 72 
Capacitors 73-75 
Capacitor 119 
Capacitor 120 
Capacitor 121 
Resistors 49-51 
Resistor 52 

Reference Designation 
Resistors 63-66 
Resistor 67, 100 and 114 
Resistor 76 

Resistors 77, 79, 31 and 110 
Resistors 78, 80 and 82 
Resistors 33-85 
Resistors 108 and 113 
Resistor 109 
Resistors 111 and 112 
Resistor 116 
Resistor 117 
Resistor 118 



1 . A method for identifying complex lithologies 
of a rock formation comprising the steps of: 



Description 

DG 502 (Tektronix) 

ICL 8038 (G.E. Intersil) 

741 (G.E. Intersil) 

10-3000 Hz 

OP 15 (Precision Monolithics) 
339 (National Semiconductor) 
TL 084 (Texas Instruments) 
IN 914 (Texas Instruments) 
White Noise (Standard 
Reference Labs, Inc.) 
0.33 uf (Kemet) 
5 pf (Kemet) 
330 pf (Kemet) 
3.3 uf (Kemet) 
0.1 uf (Kemet) 
4.99 K (Omite) 
49.9 K (Omite) 

Description 
20.5 K (Omite) 
5K (Omite) 
100 K (Omite) 
10 K (Omite) 
1.1 K (Omite) 
100 (Omite) 
320 K (Omite) 
5.1 K (Omite) 
1 K (Omite) 
25 K (Omite) 
91 K (Omite) 
15 K (Omite) 

(a) exciting at least a portion of the rock 
formation with a multi-frequency electric current; 

(b) measuring the resulting voltage in that 
portion of the rock formation; 

(c) determining the phase and amplitude of 
the measured voltage at each of the frequencies 
included in the multi-frequency electric current; 
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(d) determining the resistivity of that portion 
of the rock formation at each of the frequencies 
from the amplitude of the multi-frequency current 
at each of the frequencies and the determinations 
of phase and amplitude for the measured voltage; 
and 

(e) recording the resistivity and the phase as 
a function of frequency. 

2. A method according to Claim 1, further in- 
cluding the step of characterizing the rock forma- 
tion 

(i) as a clean rock when the resistance re- 
mains constant and the phase is negligible as a 
function of exciting current frequency; 

(ii) as a low resistivity, high surface area, 
clay-bearing rock when the resistance remains con- 
stant and the phase is a negative value which 
decreases in magnitude as the exciting current 
frequency increases; or 



(iii) a low resistivity, pyrite-bearing rock when 
the resistance decreases as the exciting current 
frequency increases and the phase is a negative 
value which increases in magnitude as the exciting 
s current frequency increases. 

3. A method according to Claim 1 or Claim 2, 
wherein the excitation is by means of an electric 
current comprised of white noise. 

4. A method according to Claim 1 or Claim 2, 
w wherein the excitation is by means of an electric 

current comprised of a plurality of discrete sine 
waves. 

5. A method according to Claim 3 or Claim 4, 
wherein all frequencies are of the same amplitude. 

rs 6. A method according to Claim 4 or Claim 5, 

wherein the plurality of discrete sine waves are 
simultaneously applied to the portion of the rock 
formation. 

7. A method according to Claim 6, wherein the 
20 stimulating electric current, l(t), is represented by 
the expression: 



I(t) = i Q Aj sin (Jt + 0.), where 

J= Ul 



lo = D.C. current component, 
A ■ gain x amplitudes, 

0 - phase-angle. 

8. A method according to Claim 4, wherein the 
excitation is by means of an electric current which 
is frequency modulated with the plurality of dis- 
crete sine waves. 

9. A method according to Claims 1 to 8, 
wherein the multi-frequency electric current com- 
prises a band width of from 0.001 hertz to 20 
kilohertz. 

10. A method according to Claim 9, wherein 
the multi-frequency electric current comprises a 
plurality of sine waves at discrete decade frequen- 
cies. 

11. A method according to Claim 10, wherein 
the discrete decade frequencies include 1 hertz, 10 
hertz, 100 hertz, 1 kilohertz, and 10 kilohertz. 

12. A method for producing a multi-frequency 
electric log of a subsurface rock formation sur- 
rounding a borehole, comprising the steps of: 

(a) traversing the borehole with a logging 
tool housing a multi-frequency electric current gen- 
erator; 



(b) introducing a multi-frequency electric cur- 
rent from the generator into the rock formation 
through a pair of in situ current electrodes as the 
tool traverses the borehole; 

(c) measuring the resulting voltage in the 
rock formation through a pair of in srju voltage 
electrodes, the measured voltage representing the 
reactive component of the electrical impedance of 
the rock formation; 

(d) determining the phase and amplitude of 
the measured voltage at each of the frequencies 
included in the multi-frequency electric current, 

(e) determining the resistivity of the rock 
formation at each of the frequencies included in the 
multi-frequency electric current; and 

(f) plotting the resistivity and phase as a 
function of frequency. 

13. A method according to Claim 12, wherein 
the logging speed of the tool is selected so that the 
low frequency end of the band width of the multi- 
frequency electric current is established in accor- 
dance with the expression: 

f = LS/60 x A p, where 



15 

f = frequency in hertz, 
LS = logging speed in m per minutes, and 
A p = depth increment in m. 

70 
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